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Abstract. We prove a global error estimate for a Lie- Trotter splitting opera- 
tor associated to the Schrodinger-Poisson equation in the semiclassical regime, 
when the WKB approximation is valid. So long as the solution to a com- 
pressible Eulcr-Poisson equation is smooth, the error between the numerical 
solution and the exact solution is controlled in Sobolev spaces, in a suitable 
phase/amplitude representation. As a corollary, wc infer the numerical conver- 
gence of the quadratic observablcs with a time step independent of the Planck 
constant. A similar result is established for the nonlinear Schrodingcr equation 
in the weakly nonlinear regime. 



1. Introduction 
We consider the nonlinear Schrodingcr equation, for t 0, 

(1.1) ied t u s + yAw £ = e a f(\u s \ 2 )u 6 . 

The function u e = u s (t,x) is complex- valued, and the space variable x belongs to 
R d . The presence of the parameter e is motivated by the semi- classical limit, e — > 0. 
Physically, s corresponds to a small ratio between microscopic and macroscopic 
quantities, so the reality is expected to be correctly described by the limit e — > 0; 
see e.g. [H] and references therein. The parameter a ^ measures the strength of 
nonlinear interactions: in the WKB regime, which is recalled below, the nonlinearity 
is negligible if a > 1, it has a leading order (moderate) influence if a = 1 (weakly 
nonlinear regime), and its influence is very strong in the regime e — > if a = 0. 
The case < a < 1 is not considered here, but it should be considered as similar 
to the case a = ([ZD- 

In this paper, we consider mostly two families of nonlinearity: 

• Nonlocal nonlinearity in the case a = 0: f(p) = K*p. 

• Local or nonlocal nonlinearity in the case a ^ 1. 

The first case includes the Schrodinger-Poisson system in space dimension d ^ 3 
(f{p) = AA~V, A e R, hence K(x) = \cd/\x\ d ~ 2 ). The second case includes the 
cubic nonlinearity (focusing or defocusing). We will also discuss why the case of 
the cubic nonlinearity is not treated in the regime a = (see Remark [ 

The initial data that we consider are of WKB type: 

(1.2) u £ (0,x) =a {x)e^°^/ e . 
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An important well-known property of this framework is related to the following 
quantities (quadratic observables), 

Position density: p £ (t,x) = \u e (t,x)\ 2 . 

Current density: J £ (t,x) = elm(u £ (t,x)Vu e (t,x)) . 

Consider the case of Schrodinger-Poisson system in dimension d 3, with a = 0. 
Formally, p e and J 6 converge to the solution of the compressible Euler-Poisson 
equation 

{d t p + div J = 0; p\ t=0 = |a | 2 , 
d t J + div (^^j + PVP = 0; J, t=0 = |a o | 2 V0 o , 
AP = \p, P(t,x), VP(t,x) -> as \x\ -> oo. 
See e.g. [6j [27] for a rigorous statement of this result. 

1.1. Time splitting methods. When simulating numerically (jl.ip . the size of e 
becomes an important parameter: if the nonlinearity / is replaced by an external 
potential V(x) (independent of u £ ), then it was proved in [24] that finite difference 
approximation requires to consider a time step At = o(e) in order to recover the 
above quadratic observables. In [3J, it was proved that these quadratic observables 
can be accurately recovered for time steps independent of e, if time splitting meth- 
ods are considered and V is bounded as well as all its derivatives. Moreover, if 
At = o(e), then the wave function u e itself is well approximated; see also [12] The- 
orem 2]. In the appendix, we extend this result to the case of unbounded potentials, 
which grow at most quadratically in space. 

In the nonlinear framework (|l.ip . numerical experiences in [3] suggest that con- 
sidering At = 0(e) is enough to recover the correct observables for time splitting- 
spectral methods, when a = 1, or a = with a defocusing nonlinearity. In the 
references mentioned so far, space discretization is considered too: in the present 
paper, we shall discuss only the time discretization, hence the above restrictions. 
In the recent paper |13j . some precise local error estimates have been established, 
showing that the assumption At — 0(e) is a sensible assumption for the local 
error to behave properly. We underscore that if crucial, the local error estimate 
is not sufficient to obtain a global error estimate, unlike in [23], because of rapid 
oscillations. 

We now briefly recall what time splitting methods consist in, in the context 
of (|1.1[) . The remark is that if the Laplacian or the nonlinearity is discarded in 
(|l.ip . then the equation becomes explicitly solvable. We denote by X\ the map 
v e (0, ■) i-> v £ (t, •), where 

e 2 

(1.4) ied t v £ + —Av e = 0. 

The above equation is solved explicitly by using the Fourier transform (defined in 
Assumption 11.21 below) . since it becomes an ordinary differential equation 

(1.5) l£ ^-^|e| 2 ^ = 0, 
hence 
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If we now denote by Y* the map u> e (0, •) i-> w £ (t, •), where 



(1.6) 



ied-, 



'tw 



= s a .f{\i 



then we remark that since / is real-valued, the modulus of w £ does not depend on 
time, hence 

(1.7) Y*w{x) = w(x)e- t£a i f ^ w ^ 2) . 

At this stage, it is already clear that whether a 1 or a < 1, the estimates for 
Y* will be rather different. We shall denote by S 1 * the nonlinear flow associated to 
CEO): 5*u £ (0,-) =w £ (V). 

We consider the Lie-type splitting operator 

(1.8) Z\ = YlX\, 

for which calculations will be less involved than for the Strang-type splitting oper- 
ator 

Zl >s = X t J 2 Y t e X t J' 1 . 
Since both X\ and Y* are unitary on L 2 , so is Z\: 
(1.9) 



-^eI|l 2 ^L 2 = I IK* ||l 2 -»-L 2 = \\Z s \\l 2 ^L 2 = 1- 



The action of Z\ on Sobolev spaces is more involved, because of the nonlinear op- 
erator Y* (in the case a = 0). In the case e = 1 with f{y) = y (cubic nonlinearity), 
the convergence of the approximate solution generated by the splitting operator as 
the time step goes to zero has been established in [3] for x £ R d , d ^ 2, and in [33] 
for x £ R 3 . 

Theorem 1.1 (From (5j[23j). Let e = 1, f(y) = y, and d ^ 2. For all uo = wf t=0 £ 

H 2 (R d ) and all T > 0, there exist C and ho such that for all At £ (0, ho], for all 
n £ N such that nAt £ [0,T] , 



(ztTu a -s 



iAt 



u 



<C(m 2 ,T) At, 



where, for j £ N, 



o ^^||zi(t)||^ (Ri) . 



Ifd = 3 and u £ H A (R d ), then 

(zf}) n u Q -S nAt u 



L 2 



<C(m 4 ,T) (At) 2 . 



Note however that these results do not directly yield interesting results in the 
case of Ijl.ljl in the semi-classical limit: in the presence of rapid oscillations as 
in (|1.1[) . the quantity rrij behaves like E~i , so the bounds in [5] [23] cease to be 
interesting. 

On a more technical level, note that even though Z\ is unitary on L 2 , the stan- 
dard Lady Windermere's fan argument, which consists in writing 

n-1 

, - u(t n ) = £ ({Z?) n - j - 1 Z^S^u - {Z? «)' 



(1.10) 



rAt^-J-l S At S jAt u ^ 
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cannot be used directly, since Z\ is not a linear operator. Therefore, nonlinear 
estimates are needed. In the case of the Schrodinger-Poisson system, the proof in 
[2"5] uses for instance the estimate 

||A _1 (uu)w;|| i 2 (R 3 ) ^ C||u|| ff l (R 3 ) ||v|| i 2 (R 3 ) ||w|| i 2 (R 3 ) . 

In the present framework, functions arc e-oscillatory (see Remark 11.41 below) , so 
the natural adaptation of the above estimates is of the form 

where C is independent of e and 

IKHin(R3) = sup (|ju e |j L 2 + ||eVu £ || L 2) 

0<£<1 

is expected to be bounded uniformly in e, unlike the standard .ff^-norm. We then 
face an e -1 / 2 singular factor in the above estimate, which ruins the approach of [23] 
in the semi-classical limit. Such phenomena explain why there is a gap between 
the proof in the semi-classical regime for the linear Schrodinger equation |12j and 
adapting the arguments of [23| to the semi-classical regime, even with the local 
error estimate of j!3| . 

1.2. WKB analysis. Given (|l.lj) with initial datum (jl.2[) . WKB method consists 
in seeking 

u e (t, x) = a £ {t, x)e^ {t ' x)/s 1 with a e « a + ea {1) + ... 
Plugging this ansatz into and ordering the powers of e, we find formally: 

O(e ): d t <P+ l -\V 



Ois 1 ): <9 t a + V0-Va 



We see that if a > 1, then the nonlincarity does not affect the pair (a, </>), which 
describes the behavior of u £ at leading order. On the other hand, if a = 1, the 
transport equation for a is nonlinear, while the equation for <f> is the same as in the 
linear case: one speaks of weakly nonlinear regime. Finally, in the case a = 0, the 
system of equations shows a strong coupling between all the terms, and is actually 
not even closed. 

In the rest of this subsection, we focus our attention on the case a = 0. An 
important remark consists in noticing that the transport equation 

(1.11) d t a + V0 • Va + ^aAcj) = -2if (\a\ 2 ) a Re (aa (1) ) , 

if it cannot be solved when is unknown, enjoys the following property: it is of 
the form D t a — ia x R, where D t stands for the vector field dt + V<f> ■ V + \A<j). 
Therefore, £>t|a| 2 = 0, and if we set (v,p) = (V<fi, \a\ 2 ), then the system in (<j>, a) 
becomes the closed system 

, . ( d t p + div(pv) = 0; p\ t=0 = |a | 2 , 

i d t v + v-Vv + Vf(p)=0; U| t=0 = V^ . 
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Note also that if we set J = pv, then (p, J) solves (|1.3[) : we have written (|1.3[) in a 
different form, which is also encountered in fluids mechanics. As a matter of fact, 
in the case of a nonlocal nonlinearity f(p) = K * p, p. lip is not correct, but since 
this term has disappeared in (|1.12l) . we do not write the correct version of (jl.lll) . 
which is a bit involved to present. In the case of a nonlocal nonlinearity, we will 
make the following assumption. 

Assumption 1.2. The nonlinearity f is of the form f(p) = K*p, where the kernel 
K is such that its Fourier transform, defined by 

satisfies: 

• Ifd^2, 

sup (1 + |e| 2 )|i?(OI < oo. 
£eR d 

• Ifd^ 3, 

sup |£| 2 |£(£)| <oo. 

£GR d 

Typically, this includes the case of Schrodingcr-Poisson system if d ^ 3, where 
f(p) is given by the Poisson equation 

Af = \p, f, V/ -> as |x| -> oo, 

with A £ R. This equation can be solved by Fourier analysis if d ^ 3 (K(£) — 
— A|£| -2 ); if d ^ 2, this is no longer the case, as discussed in [25]. Under this as- 
sumption, (|1.12|) has a unique solution (v, p) £ C([0, T];H S+1 x (H s nL 1 )) provided 
that the initial data are sufficiently smooth, with s > d/2 + 1, from [15] (see also 
[HUH [57], an d Section [3] for the main steps of the proof). 

Proposition 1.3. Suppose that f satisfies Assumvtion \1.2l Let ao,<f>o £ S'(R d ) 
with (V(f>o,ao) £ H s+1 x H s for some s > d/2. There exists a unique maximal 
solution (v,p) £ C ([0, T max ); H s+1 x (H s n L 1 )) to (pTT2]) . In addition, T max is 
independent of s > d/2 + 1 and 

T max <+oo=^ / (||«(t)|| W ri.- + ||a(*)|| W 'i.oo)tK = +oo. 
Jo 

Remark 1.4. We note from (|1.12p that even if no rapid oscillation is present initially 
in (]1.2I) . then ui t= o = and 9t^|t=o 7^ 0, so the solution u £ is not e-oscillatory at 
time t = 0, but becomes instantaneously e-oscillatory. 

We emphasize the fact that under Assumption 11.21 and for fixed e > 0, given 
Uq £ L 2 (R d ), (|1.1[) has a unique, global solution u e £ C([0, 00); Li 2 ). Moreover, 
higher Sobolev regularity is propagated globally in time (the nonlinearity is L 2 - 
subcritical) ; see e.g. [10] . 

1.3. Main results. Our main result measures the accuracy of the time splitting 
operator so long as the solution to (|1.12[) remains smooth. 

Theorem 1.5. Suppose that d 1, a = in (jl.ll) . and that f satisfies As- 
sumption FTM Let {<j>o,a ) £ L°°(R d ) x H s (R d ) with s > d/2 + 2, and such 
that V0o £ H s+1 (R d ). Let T > be such that the solution to (|1.12j) satisfies 
(v,p) £ C([0,T];H s+1 x H s ). Consider u e = S^uf, solution to ^B) and Uq given 
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by (|1.2j) . There exist eo > and C, Co independent of e € (0,£o] such that for all 
At G (0,co], for all n G N ,suc/i t/iat t„ = nAt G [0, T], the following holds: 

1. There exist <f> £ and a e with 

sup (||a e (t)|| HS(Rd) + ||V0 £ (t)|| ffs+ i (Rd) + U e {t)\\ L ^ (Rd) ) < C, Ve e (0,£ ], 
te[o,T] 

sucfc t/iat u e (t,:z;) = a e (t, a;) e ^ E (M)/s for all (t,x) G [0,T] x R d . 

2. There exist <f> £ n and a e n with 

IKIIff»(R*) + l|V^|| ff .+i(Rd) + HdU-(Rd) s$ C Ve € (0,£ ], 

suc/i that (Z^*)™ (aoe l ^°/ e ) = a^e l ^"/ e ; and the following error estimate holds: 

\K - o'COIIff.-! + l|V0= - V<f (t„)|| ffs + ||0» - <j> e (t n )\\ LOO < CAt. 

Note that in the above result, the phase/amplitude representation of the exact 
solution u e and the numerical solution is not unique. We infer the convergence of 
the wave functions in L , by reconstructing the numerical wave function: 



Corollary 1.6. Under the assumptions of Theorem I J. 51 there exist Eq > and 
C, Co independent of e G (0,£o] suc/i that for all At G (0,co], /or all n € N such 
that nAt G [0,T], 

|| (Zg *) n UQ — 5*"Uo|| i2 ^ R ^ < 
We also get the convergence of the main quadratic observables: 

Corollary 1.7. Under the assumptions of Theorem \1.5l there exist e > and 
C, Co independent of e G (0,£o] suc/i that for all At G (0,Co] } /or all n £ N such 
that nAt G [0,T], 



|(Z £ Ai )Xr-|p £ (*«)| 
Imf £ (Z^)™ U §V(Z e At )X) 



2 



< CAt, 

L 1 (R d )ni°°(R d ) 



< CAt. 

L 1 (R d )nL°° (R d ) 

These results may seem limited, inasmuch as they address only a specific regime, 
and say nothing on the large time behavior. We emphasize the fact that the behav- 
ior of u e as £ — > at time where the solution to (|1.12j) ceases to be smooth is still 
an open question. Therefore, the analytical tools to analyze the splitting operators 
are missing, due to a lack of precise estimates on the exact solution. Typically, all 
the results presented here highly rely on the fact that a WKB regime is considered. 

1.4. Weakly nonlinear regime. We now consider the case a ^ 1 in ([l.ljl . which 
turns out to be quite easier to treat. To begin with, the assumption on the nonlin- 
earity is weaker, and we allow local interactions. 

Assumption 1.8. The nonlinearity f is of the form / = /1 + /2, where f\ satisfies 
AssumptionEM and f 2 G C°°([0, 00); R+) 7 with / 2 (0) = 0. 

Remark 1.9. The assumption /2(C)) = is here merely to simplify the presentation, 
since replacing / with / — /2(C)) in (jl.ip amounts to replacing u E with u E e lt ^ 2 ^-°^ e . 

Proposition 1.10. Suppose that d ^ 1, / satisfies Assumption ! 1 . 81 and that a 1 
in (fTTT]) . Let Oo,a ) G H s+2 x H s with s > d/2 + 2. Let T > be such that the 
solution to 

9t0+i|V0| 2 =O; 0| t=o = 



TIME SPLITTING FOR SEMI-CLASSICAL NLS 



7 



satisfies <fi G C([0, T]; H s+2 ). Consider u £ = 5*Uq solution to (|1.1[) wit/i a ^ 1 arte? 
Mg given by (|1.2|) . There exist eq > and C, Co independent of £ € (0,£o] swc/i £/ia£ 
/or aZ/ At G (0, Co], /or a!/ n £ N smc/i i/iai nAt G [0, T], i/ie following holds: 

1. //we set a e = ife' 1 ^^ , then 

sup ||a e (t)||H S(R£i ) ^ C, VeG (0,e ]. 
fe[o,T] 

2. There exist 4> e n and a e n with 

IKIIho(r") + H^rJIffs+W) <C, Vee (0,e ], 
sitc/i t/iat (Z^ t ) n (aoe^°/ e ) = a e n e % ^ n l e , and the following error estimate holds: 
IK - a e (UH ff3 - 2 + HI ~ <t>{t n )\\ H s < CAt. 

In particular, 

\\{Z^) n ul-S^ At ul\\ L ^C^. 

It may seem surprising that even in the weakly nonlinear regime a = 1 , the result 
is local in time, and valid only before the possible formation of caustics. As a matter 
of fact, the behavior of the nonlinear solution u e is essentially not understood past 
the caustic; see e.g. [7]. 

Notations. Throughout the text, all the constants are independent of e G (0, 1]. 
For (a e )o<e^i and (/3 e )o<£^i two families of positive real numbers, we write a e < j3 e 
if limsupa 6 //3 e < oo. 

2. Action of the numerical scheme in the WKB regime 

Our approach consists in sticking to the WKB framework. We write the solutions 
to under the form 

(2.1) a E (t,x)e l ^^/ E , 

with a £ and <\f bounded in H s (R d ) uniformly in e G (0,1]. Here, the "phase" 
<j) e is real-valued, and the "amplitude" a e is complex-valued; of course, such a 
representation is not unique. As a matter of fact, both a e and <fi 6 must be expected 
to depend on e. Functions of this form will be referred to as WKB states throughout 
the text. From now on, and up to Section [3 we assume a = in (|l.ll) . 

2.1. A stable phase/amplitude decomposition. The important remark con- 
sists in noticing that the flows associated to (|1.5[) and (|1.6[) preserve the structure 
of WKB states. 

Nonlinear flow. In the case of the nonlinear flow (|1.6[) . the exact formula (|1.7p 
shows immediately that a WKB state evolves as a WKB state: if wf t=0 = a £ e llpe l e , 
then the solution to (|1.6[) is given by 

w%t,x)=a%x)e<^^- t ^ aE ^ 2 ))^. 
This is indeed of the form (|2.1[) . with 

a £ (t,x) = a s (x), ^(t,x)=cp £ (x)-tf(\a £ (x)\ 2 ). 



8 



R. CARLES 



We can therefore rewrite the action of Y" e * on WKB states as the action of the flow 
yf on phase/ amplitude pairs (cj>, a) characterized by 



(2.2) 




= 0; 0r t=o = (f> , 



t=o 




Linear flow. The analysis of the linear flow (|1.4|l is less straightforward, and requires 
more care than the nonlinear flow. Consider the system 



(2.3) 



Note that this is not exactly the system corresponding to standard WKB analysis, 
because of the term eAa £ in the second equation, which is discarded in WKB 
approximation. The first equation is an eikonal equation, which has a smooth 
solution at least locally in time (see e.g. [7]), and energy estimates then follow 
easily for the second equation. We emphasize the fact that (|2.3j) is equivalent to 
(|f ,4j) in the case of initial WKB states (|2.1|1 . at least locally in time, modulo the 
eikonal equation. Indeed, given an initial phase <^q, we can solve, locally in time, 
the eikonal equation 

(2.4) ^ + i|v^| a = 0; <j>\ t=0 = 4>l- 

In general, the solution to (|2.4p does not remains smooth for all time, due to the 
formation of caustics (see e.g. [7]). We note that w e = V(f> e solves a (multidimen- 
sional) Burgers equation 

9(!o e +ic e -Vro E = 0; w* t=0 = V<j>%. 

This remark will be used to derive a priori estimates for the system (|2.3p . Once (jf 
is known, then v e , solution to (|1.4[) . and a e , are related through the formula 

which yields an obvious bijective correspondence between these two functions. Even 
though there is no uniqueness, we conclude that if the initial datum is a WKB state, 
«u_ = age^o' 5 , then at least locally in time, v e remains a WKB state, since it 

can be written as v £ = a e e l ^^ e , where ((jf,a e ) is the solution to (|2.3j) . Following 
the same convention as in the case of the nonlinear flow, we denote by X\ the flow 
acting of phase/amplitude pairs, 



Xl 



'■(t) 
'■(t) 



where (4> £ , a £ ) is the solution to (|2.3|) . Similary, we write Z\ — y\X\. 

2.2. Rewriting the splitting method in the WKB regime. Instead of ana- 
lyzing directly the equations (|1.4|) - (|1.6|l . we shall work on (|2.3|) — (|2.2[) . in view of 
the previous subsection. We denote by IF the wave reconstruction operator 
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and we note the identity, which is the key conclusion of the above analysis: 
(2.5) WZt W = Z% (aV^) . 

In view of the obvious remark 

we see that working with (<jf ,a £ ) is not equivalent to working with the wave func- 
tion a £ e l ^ c / £ . However, we only use the fact that the numerical solution can be 
represented by this decomposition, and no uniqueness argument is needed, except 
the fact that the solutions to (|1.4[) and (|1.6[) , respectively, are unique. 

We finally notice that the form (|2.1[) (with a £ and <j> £ bounded in H s (Il d ) uni- 
formly in e) is preserved by the exact flow. This is so thanks to the gauge invariance 
of the nonlinearity, as noticed in |16| (an aspect which also appears when solving 
(|1.6[) ). In the case of a local defocusing nonlinearity (typically f(p) = p), recall 
the original idea of Grenier [17] to study the semi-classical limit for seek 
the solution u £ to (|1.1[) under the form (|2.ip . also with (f> £ real- valued and a £ 
complex- valued. One gains a degree of freedom, and the choice of Grenier consists 
in imposing 



(2.6) 



dtf + I|V<f | 2 +/(|a E | 2 ) = 0; 0f t=o = ^ 
1 e 

d t a £ + V0 e • Va £ + -a £ A(j) £ = i- Aa £ ; a[ t=0 = a . 



This choice differs from the standard Madelung transform, which is limited by the 
presence of vacuum (zeroes of a £ ; see [5]). Under Assumption 11.21 the adaptation 
of the approach of Grenier can be found in [TJ [5T| (see also [22] HS1 US] f° r the 
case of Schrodinger-Poisson system in low dimensions, where low frequencies cause 
technical difficulties). We denote by S £ the flow associated to (|2.6[) . 

Remark 2.1. In (|2.6|) . the initial data are supposed implicitly independent of e. 
This is merely for the sake of consistency in future references. Throughout this 
paper, the flow associated to (|2.6|) will be considered for initial data which may 
depend on e, but which are uniformly bounded in suitable Sobolev spaces. 

Instead of analyzing directly the splitting method for (|1.1[) as presented in Sec- 
tion 11.11 we shall therefore analyze a splitting method for (|2.6|) : when the term / 
is discarded, we recover (|2.3[) . which is solved alternatingly with (|2 . 2|) . The lat- 
ter system consists indeed in dropping out the Laplacian in (|2.6p . since all spatial 
derivatives have disappeared. 

3. Technical background 
As noticed in [1], the following lemma turns out to be very helpful. 
Lemma 3.1. Let s 0. Under Assumvtion \1.2[ there exists C such that 

(3.1) \\vf( P )\\ H , + . {m < c (IMU*(r«) + IMUW , v P G H s (R d ) n l 1 ^). 

If in addition s > d/2, there exists C such that 

(3.2) \\f(p)\\ L ^) < C (|M|^ (R «) + \\p\\l^)) , Vp G H s (R d ) n i 1 (R d ). 
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Proof. By Plancherel formula, 



v/(p)ii^ +1(Rd) = / \z\ 2 {i+\ey +1 \K(t)\ 2 \m\ 2 dt 



^ sup (l + |£| 2 )|A^)| \\ P \\h, 
hence (a weaker version of) the lemma in the first case of Assumption 1 1.21 If d ^ 3 



Ki 2 (i + ici 2 ) s+i \m\ 2 \p(o\ 2 dt < sup \t\ 2 \Km / ier 2 ip(6i 2 ^ 

l?Ki \e£R d / -/ISKi 

< C\\p\\l-w [ r d - 3 dr^C\\p\\l 1{Rd) , 
Jo 

where we have used spherical coordinates and Hausdorff- Young inequality. This 
yields the first part of the lemma. For the second part, we use the same tools, 

||/(p)|U~ < (2^-^11/^)11^ = \\Kp\\ L ,. 
Split the integral between the two regions {|£| ^ 1} and {|£| > 1}: 

\m)\\p(0\dt<C\\p\\ L ~ [ r d - x K(r) dr<\\p\\ L1 , 

Jo 

\K(0\\mm < < since 3 > d/2. 

'I?l>i 

This estimate is not sharp, since we do not use the decay of K at infinity. □ 

As in PQ, we infer the following result, concerning the exact solution, that is, the 
solution to (|2.6p . This result implies the first point of Theorem 11.51 



Proposition 3.2. Suppose that d 1, and that f .satisfies Assumption ll . e A Let 
(V<fio,ao) G H s+1 x H s with s > d/2 + 1, and let T > be such that the solution 
to (fTTT2|) satisfies (v,p) G C([0, T]\H S+1 x H s ). Then there exists e > such 
that for all e G (0,eo\, (|2.6[) has a unique solution, which satisfies (V</> e ,a e ) G 
C([0,T];H s+1 xH s ), uniformly in e G (0,e ]: ^ere exisis C(T, ||a ||/f», ||V0 o ||jj-+O 
independent of e G (0,£o] such that 

sup (\\a £ {t)\\ H s (Rd) + \\V(f> e (t)\\ H s + i {Rd) ) s$ C (T, Haolli^Rd), || V^oll^+^Rd)) . 
*e[o,T] 

If in addition O € L°°(R d ), then cj) E G C([0, T]; L°°(R d )) and 

sup ||0 e (i)|Uo O(R .< 1) < ||0o||l~ + G(T, ||a || ffs(R£i) , ||V0 o ||ff=+i(Ri)) ■ 
te[o,T] 

Sketch of the proof. Let w 6 = V0 e . By differentiating in space the first equation 
in (|2.6p . we see that any solution to (|2.6p must solve 

C d t w £ + w e • V» E + V/ (K| 2 ) = 0; u;f t=0 = V0 O , 

(3-3) { i e 

y d t a £ + w e ■ Va £ + -a £ div w e = i-Aa £ ; a jt=o = a o- 

The left hand side corresponds to a hyperbolic symmetric system for the unknown 
(w e , Rea e , Ima e ) G R d+2 , thanks to Lemma l3~T| and the shift in regularity between 
w £ G iP +1 and a e G IP. The right hand side of (|3.3|) is a skew-symmetric term, 
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which does not appear in H s energy estimates. The key point to notice is that unlike 
what would happen in the case of the nonlinear Schrodinger equation, the terms 
V/(|a e | 2 ) and a £ drvw e are not quasilinear, but semilinear (they can be treated as 
perturbations) , in view of Lemma 13.11 and the functional framework. By standard 
theory (see e.g. [2]), (|3.3j) has a unique solution (w e , a £ ) £ C([0, r]; H s+1 x H s ), for 
some r > independent of e G (0,1]. 

We can take r ^ T for e sufficiently small. Indeed, if T' denotes the lifespan 
of (|3.3|) in the case e = 0, then necessarily T' > T, for if we had T" ^ T, then 
by uniqueness for the Euler-Poisson system, \a\ 2 = p e C([0,T]; H s n L 1 ) and 
«; = y e C([0,T]; iT s+1 ). Back to the transport equation in (|3.3[) . we infer that 
a e L°°([0, T]; which yields a contradiction. 

Finally, we note that a e , u> e and e are related through the formula 
d t <j> £ + l\w°\ 2 + f(\a s \ 2 ) =0; 0f f=o = 0o 

Therefore, if u> £ and a e are known, then <p e is obtained by a simple integration in 
time, and the last estimate of the proposition follows from (|3.2|) . □ 

Note that the above result is expected to be valid only locally in time, since the 
solution to (|1.12l) may develop a singularity in finite time. In that case for fixed 
e > 0, a £ may become singular, or stay smooth but become e-oscillatory for large 
time, as suggested by the simulations in [5]. This can be understood as follows: 
for large time, several oscillations are expected in u E , so they cannot be carried by 
only one exponential function as in (|2.1[) . therefore, a £ is rapidly oscillatory, and 
its £P-norm is not bounded uniformly in e £ (0, 1]. 

The analysis of [T] also implies the following result. 

Proposition 3.3. Let d 1, and f satisfying Assumption 1 1.8l Let R > and 

s > d/2 + 2. There exists T = T(R) > such that if 

\\a>o\\H°(R d ) + l|V0o||_H-s+i(Rd) < R, 

then (|1.12l) has a unique solution (v,p) £ C([0, T]; H s+1 x H s ). There exist £q > 

and K = K{R) independent of e € (0,£o] such that if in addition (\7tpo,ba) G 
H s+i x H s satis fi es 

t/ien /or a?Z t g [0, T], t/ie solutions to (|2.6p initial data (4>o,ao) and (tpa,bo), 
respectively, satisfy: 

\\a E (t)-¥(t)\\ Hs +\\V0 E (t)-Vip e (t)\\ H s+i < A"(||a - 6 || ffs + ||V0 O - V^o|| H .+i) . 
There exists k = k(R) such that if in addition (j>o, ipo £ L°°(R d ), then 

\\<j) £ {t) - Lp e (t)\\ L ^ < H0O - ipolk- + « (||oo - &o||h« + || V0 O - V^o||^+i) • 

4. Estimating the approximate flow 

In this section, we prove various estimates concerning the flows involved in the 
definition of the numerical scheme, X\ and y\ . 
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4.1. The nonlinear operator. Unlike what happens in most cases when studying 
splitting operators, the most delicate operator to control is the linear one, denoted 
here by X*, while in the present framework, y l e turns out to be the simpler of the 
two. 

Lemma 4.1. Let s > d/2 and <j> £ ,a £ G S'(R d ), with (V$j,a§) G H s+1 x H s , for 
some s > d/2. The solution to (|2.2|) is given by 

r(*) = 0S-*/(l4| 2 ); a £ (t)=af ) . 

In particular, there exists C = C(/i) such that if ||aol!i°° ^ M> 

IKWII*' = |K|| H s, ||V0 E (t)|U« +1 < ||V^||^ +1 +Ct||^|| ff3 , Vt^Q. 

Finally, if 4>q G L°°(R d ), iften i/iere exists C = C(fi) such that if || a§ || ^ fJ>, 

\\4> e (t)\\L°° sS \\<t>o\\L°° + Ct\\al\\ H s, Vt ^ 0. 

Proof. Since s > d/2, _ff s (R d ) is a Banach algebra embedded into C(R d ), hence 
the formula for <p e . The estimates are straightforward consequences of Lemma |3.1[ 
and of the tame estimate jj» < HsII-F/ 3 + II /II if* IMIl 00 ■ O 

4.2. The linear operator. Wc now consider (|2.3j) . The following lemma is a 
variant of [HI Lemma 3.2]. 

Lemma 4.2. Let s > d/2 + 1 and \i > 0. There exists t = r(/i) > such that if 

||vo||j?» sS Mi 

i/ien i/ie (multi- dimensional) Burgers equation 

(4.1) ftv + v ■ Vv = 0; V| t=0 = v 

ftas a unique solution v G C([0, r]; which satisfies 

sup ||v(*)||j?-. s$ 2/i. 
te[o.r] 

Proof. Local existence of a unique H s solution follows from a global inversion the- 
orem (see e.g. [7]), so we focus on the energy estimate. We have 

\^Mh° = (v,a t v> H3 - (A*v,A s 9 t v) i2 = -(A s v,A s (v • Vv)) L2 , 

where A = (1 — A) 1 / 2 . Introduce the commutator 

\j t M 2 H° = v • VASv >l 2 + (A s v, v • VA s v - A s (v • Vv)) L2 . 

By integration by parts, the first term is controlled by 

|(A s v,v VA'V) L2 | ^ illvllir.ll divv^- < 

where we have used Sobolcv embedding and the assumption s > d/2 + 1. The last 
term is estimated thanks to Kato-Ponce estimate [20] 

(4.2) \\A°(fg) - fA s g\\ L * < \\Vf\\L-\\g\\ H .-i + ll/IMMU-, 
with / = v and g = Vv: 

|(A s v,vVA s v-A s (vVv)) i2 | ||v|| H . ||A S (v Vv) - v VA s v|| i2 

<l|Vv|| LO c||v||^ <||v|&- 
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We infer Jj||v||fjs ^ C|| v lHf=: an d the result follows by comparing with the ordinary 
differential equation y = Cy 2 . □ 

Lemma 4.3. Let s > d/2 + 1 and \i > 0. If the solution v to (|4.1[) satisfies 

||Vv(t)||ioo OsCtsCr, 

then there exists c independent of (i and r such that 

sup ||v(i)|| ffs < e c ' lt i|v(0)|| ffs , < t < r. 
te[o,r] 

Proof. This lemma is a straightforward consequence of the tame estimates used in 
the proof of Lemma 14.21 □ 



Proposition 4.4. Let a > d/2, fx>0. Suppose that (V^ag) e H cr+1 x H cr , with 

IIV^oIIj^+i 5$ /i, \\a e \\ H ° < 

There exists r = r(/i) independent of e such that (|2.3p /ias a unique solution, with 
(V0 £ ,a e ) S C([0,t];F <7+1 x and 

sup ||V0 £ (i)|| H „+i «S 2/i, sup ||a e (i)|| H . < 2//. 

te[0,r] i£[0,r] 

If in addition <pl £ L°°(R d ), then <p € e C([0,t];L°°) and 

sup ||<^(t)|| L ^ ^ ||0olU^ + r A«- 

te[0.r] 

Proof. From Lemma l4~2l (|4.1[) has a unique solution v <G C([0, r]; H a+1 ), such that 
V| t=0 = V0g, with ||v(t)|| ff „+i ^ 2^i for t G [0, r]. Now let 

e (*) = <^o - \f \v(<r)\ 2 d*. 

We note that <9t(v— V(/> e ) = dt~v—Vdt(jf — 0, so v = V0 e , and the result concerning 
£ follows. 

The existence of a solution a £ follows for instance from the fact that it is given 
by a E = v e e~' l< ^' l e , where v 6 € C(R;H a ) is the solution to the linear Schrodingcr 
equation (|1.4p with initial datum a§e l< ^°' e € H a . So we are left with the energy 
estimate: since iA is skew-symmetric, 

~hm, = (A^,(d t -i £ -A 



V<f • Va £ + ^a e A(jf 



L 2 



By integration by parts, 



A CT a e , V<f • VA"a e + — A a a £ Acf> E \ = 0, 

2 / L 2 

so we have 

+ i (A CT a £ , A a a 6 A<fi £ - A a (a E A0 £ )) L2 . 
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Kato-Ponce estimate (|4.2|) for the first line, and tame estimates for the second line 
then yield 

J t W ae \\H° < ll« £ ||^(!|a £ ||^l|V 2 ^||L» + ||V ( /» £ ||^||Va||L- 
(4 ' 3) + \\A^\\ L ~\\a £ \\ H ° + ||A0 8 ||^||o||l-) 

< \\v<t> e \\ H « + 4a e \\ 2 H „, 

since a > d/2, and the result follows from Gronwall lemma, by decreasing r if 
necessary. □ 

Remark 4.5. The above proof suggests that the shift in regularity, between <j) E 
and a e , cannot be avoided. Note that this phenomenon shows up when the free 
Schrodinger equation (|1.4[) is solved, in terms of WKB states. 

Proposition 4.6. Let a > d/2, fi > 0. Suppose that the solution to (|2.3p satisfies 

\\V(f> e {t)\\ W i,ao \\a £ (t)\\ w i. x ^ O^t^T. 

There exists c independent of e, fi, r such that the solution to (|2 .3[) satisfies 

||V<f + \\a'(t)\\ H , < e c ^ (||V0§|| H - + i + ||o§|| H , O^Kr. 

Proof. Lemma 14.31 readily implies 

||V<^(i)|| H „+i e^HV^H^+i, t < r, 



for some c independent of e, /x, r. Back to the proof of Proposition 14.41 simply 
apply Gronwall lemma to (|4.3p . □ 

4.3. The splitting operator. In view of Lemma 14.11 and Proposition 14.41 we 
readily have: 

Corollary 4.7. Let s > d/2, fx > 0. Suppose that (V^a 5 ) e H s+1 x H s , with 

l|V$)llir*+ I < IKII-Er 5 < M- 
There exists r = > independent of e such that Z\ \% \ = l ' V with 



sup ||V0f lljja+i ^ sup ||of < 4/x. 

tG[0,r] tS[0,r] 

// m addition (jf € L°° , with \\4>q\\l°° ^ M> ^en, itp to decreasing r, we have 

SUp ||0|||ioo 4yU. 
t£[0,r] 

Proo/. Lemma IO implies that y\ (^j = with 

llotllff- = lloollff'. l|V^|| ff a+i < fi + Ct, t ^O. 

We then apply Proposition ^. 41 with cr = s. We note that the L°° regularity for the 
phase is propagated by both operators, and the estimate follows easily. □ 

In view of Lemma l4.1l and Proposition 14.61 we also infer: 
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Corollary 4.8. Let s > d/2, fi > 0. Suppose that Z\ f^°J = (tlj> with 

||V0f || w i,oo < ||o| < M. O^i^r. 

T/iew iftere exists c independent of e, fi, t, such that 

||V#|| H « +1 + s; e c ^ (||V0g|U.+i + K||h0 , ^ t < r. 

5. Local error estimate 

We recall the result (and resume the notations) from [T3] concerning the local 
error estimate in the context of For an operator A, possibly nonlinear, wc 

denote by £ A the associated flow: 

d t £ A {t,v)=A(£ A (t,v)); £ A (0,v) = v, 

The results presented in this section rely heavily on the following result. 

Theorem 5.1 (Theorem 1 from PS])- Suppose that F(u) = A{u) + B(u), and 
denote by 

5*(u) = £ F (t,u) and Z\u) = £ B (t,£ A (t,u)) 
the exact flow and the Lie- Trotter flow, respectively. Let C(t,u) = Z (u) — S (u). 
We have the exact formula 

C(t,u)= f f * d 2 £ F (t-T 1 ,Z T i(u))d 2 £ B (T 1 -T 2 ,£ A (T 1 ,u)) 
Jo Jo 

x [B, A] (£ B (t 2 ,£ a (n,it))) dr 2 dn. 

We emphasize the fact that in [T3] , this result is established for general operators 
A and B. In particular, both operators may be nonlinear. In the case of (|1.4p - (|1.6[) . 

A = i £ -A; B(v) = --f(\v\ 2 )v; F(v) = A(v) + B(v). 

We have omitted the dependence upon e in the notations for the sake of brevity. 
The linearized flow 8 2 £f is characterized by d 2 £F(t,u)uiQ — w, where 

'- + — Aw = / (|u| 2 ) w + f (uw + uw) u; w\t=o = wq. 

We note that it is not compatible with our approach, inasmuch as it does not 
preserve the (monokinetic) WKB structure: if u = ae l ^l £ and wo = boe lipa / £ , then 
the equation becomes 

iedtw + yA» = / (H 2 ) w + f {ae- l4,/e w + ae #/e w) ae^ /e ; w {t = = b e tVo/£ . 
In general, this is not compatible with a solution of the form 

w = b £ e ilfi * /s , 

with b £ and tp £ uniformly bounded in Sobolev spaces. Possibly, w should rather be 
seeked as a superposition of WKB states, 

3 

Another, less technical, way to see that the local error should not be expected to be 
a single WKB state consists in going back to the definition. We have seen that the 
numerical solution remains of the form (at time t n = nAt) u £ n {x) = a £ n {x)e l ^ x ^ 1 ' £ , 
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while the exact solution is of the form (Proposition ^. 2j) u e (t, x) = a e (t, x)e 1 ^ ^' x '' e . 
Thus the local error is 

C{t n ,u )(x) = a^xy^l" - a s (t,xy<f s ^ £ , 

and it is very unlikely that this can be factored out as 

C(t n ,u )(x) =<(x)e^"^/ £ , 

with a e n and tp e n uniformly bounded in Sobolev spaces (consider for instance the 
trivial example, L = (e ixi ' e — l) e - '^"). 

This aspect is another motivation for working with the system (j2.3p ~ (|2.2p instead 
of the standard one (|1.4[) - (|1.6[) . We therefore consider the operators A and B 
defined by 

We note that with this approach, neither A nor B is a linear operator. 

Lemma 5.2. Let A and B defined by (|5.1[) . Their commutator is given by 

I A Ri (<i>\ - ^V0-V/(|a| 2 )-div/(|a| 2 V0)-£div/(Im(aVa)) 
[A ' Bl \aJ-\ Va.V/(H 2 ) + iaA/(|a| 2 ) 

As a consequence, if s > d/2, \\ X/<j)\\ H s+-2 ^ M, ||a||#s+i ^ M, then there exists 
C = C{AI) independent of e € (0, 1] such that 

[W*Wf|. -ak /M-»<C(|v#| HM . + w^) 1 

W V 6 / \ WHh- < C7||o|| ff .+i. 

In particular, 

\\<p\\l~ ^C(||V0|| ff3+2 + || a ||^ +1 ). 

Proof. By definition (see [13l Section 3]), 

[A,B]v = A'(v)B(v) - B'(v)A(v). 

We have, since / is linear in its argument, 

, U\ (<p\ f -V0 • \ 

W ^6/ ^-V^ ■ V6 - Vip ■ Va - ±oA0 - ±aA</j + i|A6 J ' 

fl/ ^ M = f-f{ab + ab)^ = /-2/(Re(a&)) 
We compute 

B' ^ i4 ^ = ( 2/ (Re ' Va)) + / (l a ' 2A ^) + £f ( Im ^ Aa ^ 

The main point is then to notice the factorizations 

2Re(aV(j) ■ Va) + |a| 2 A0 = div (|a| 2 V0) ; Im(aAa) = divlm (oVo) , 

and to recall djf(p) — f(djp), 1 j ^ d. 

The estimates of the lemma then follow from the explicit formula for [A, B], from 
the fact that H s+2 (R d ), H s+1 (R d ) and H s (R d ) are Banach algebras, from ([3TT]) . 
and from the embedding i/ s+2 □ 
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We have the explicit formula 

_ /V-*/(M a ) 

a 



and we readily infer 
(5.2) d 2 £ 
Finally, we compute 



(p\ _ fip - 2t Re / (ab) 



bo) \b(t)J> 

d t ip + \7(j)- V(^ + 2Re/(a6) = 0; ip\ t=0 = tp , 

(5.3) { l £ 

d t b + V0 ■ Vb + Vip ■ Va + - (bA(f> + aAtp) = i - Ab; b\ t=0 = b . 

Lemma 5.3. Let s > d/2. Assume that (V0, a) e L 1 (I;H S+2 x H s+1 ), where 
€ I. There exists C independent of e <G (0, 1] such that if (Vipo, bo) € H s+1 x H", 
the solution to (|5.3p satisfies for all t G 7, 

||&(*)ll*- + II + i SS {\\b Q \\ H . + HV^oll^ + O e C J"o(ll««llH.+i + IIV*(r)|| H . +a )«*r. 

// m addition ipo £ L°° , then 

< (IIVoIIl- + ||6o||ff. + ||VV0|| ff . + i) e C /o(ll«WII«. + i + IIV*MI| H . +a )<ir_ 

Proof. Set w = V</j: (|5.3|) implies 

<9 t u; + V0 ■ Vi« + V 2 ■ u; + 2RcV/(a&) = 0; w\t= = V</? , 

(5-4) { 1 e 

d t b + V0 • V6 + w ■ Va + - (6A0 + adiv w) = i-A6; 6| t=0 = b . 



As in the proof of Proposition ^. 2[ the term iAb being skew-symmetric, it does not 
show up in energy estimates. Using Lemma I3.R we have the estimate 

\\w(t)\\ H s +1 + \\b(t)\\ H s <; \\w \\ H s +1 + \\bo\\H- 

+C ( (||V0(T)|| ff . +a + ||o(r)|| H . + i)(||«;(T)|| ff . + i + ||6(r)|| H .)^. 



and the first estimate of the lemma stems from Gronwall lemma. 

The second estimate then follows from the first equation in (|5.3[) (integrated in 
time), and ([372]) . □ 

Putting these estimates together, and using Theorem 15.11 we obtain a result 
which is crucial in the proof of Theorem 11.51 

Theorem 5.4 (Local error estimate for WKB states). Let s > d/2 + 1 and \i > 0. 
Suppose that 

\\Vcb £ \\ HS+ i ||a e || ff3 s^. 

There exist C, cq > (depending on \i) independent of e € (0, 1] such that 

4(t>=<i)-<t>C:!;r 

where A e and \I/ e satisfy 

||Vtf e (t)||jf. + \\A*(t)\\ H .-x < Ct 2 , Q^t^c . 
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If in addition 1 1 <^> e 1 1 i °^ ^ (J-, then (up to increasing C) 

W{t)\\L~> ^Ct 2 , O^t^Co. 

Proof. Let t G [0, c], and fix ti,T2 such that ^ r 2 ^ n ^ t. Introduce the 
following intermediary notations: 



£a In, 



f! 1 



^1 



[B,A] 



*2, 

9 2 ^b - r 2 , 
£b I n 



'i II _ I n 



2 I _ I 



£a n 



Then in view of Theorem 15.11 we have 

In view of Proposition FOl we have, uniformly on [0,c], for c sufficiently small, 

||V#||jr.+i <2/i, ||o?||jr. <2/i. 
Now Lemma 14.11 implies (up to decreasing c) 

||V$j||h*+i < 3//, ||a|||H<» s% 3/z. 
From Lemma 15.21 we infer 

||V01|| H . <4fi, ||o||| fi .-i < 4/x, 
provided that s — 1 > e?/2. In view of (|5.2j) . we have 

al = <4 0I = ^-2( n -r 2 )Re/(afa|), 

and therefore 

l|V</>||| ffs -i < 5/i, ||oI||ji— i s$ 5/a, 
since s — 1 > d/2. Now Corollary 14.71 implies 

\\V4> e 2 \\ H s+i ^ 4ft, \\a s 2 \\ H s <4/x. 

Finally, Lemma 15.31 yields, up to decreasing c one last time, 

d 2 £ F (t - n, = (J) , with \\vr\\ H . < io M , H^IU-i < io M . 

The first estimate of the theorem then follows by integrating with respect to (ti, t 2 ) 
on {0 ^ r 2 ^ ri ^ £}. The L°°-estimate of the phase follows similarly. □ 

Back to the wave functions, we obtain an estimate similar to the one presented 
in [H Section 4.2.2]: 
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Corollary 5.5. Let s > d/2 + 1 and \i > 0. Let </>g <E L 00 ,^ e H s with 

||$>lk°° < Mj I|V0o||jj»+i < ft, \\a E Q \\ H - < /«. 
There exist C, cq > (depending on fi) independent of e G (0,1] sitc/i £/ia£ 



(a§e 4 «/ e ) - S£ (age**' 6 ) 



^ C— , < i < co. 

L 2 £ 



Proof. We have 5**Mq = a e \t)e 1 ^ '® ,£ where of and 4> £ are given by Proposition 
and 

a% - a £ (t) = A £ (t), # - = ¥«(t), 
where A £ and \I/ e are given by Theorem 15.41 We compute, since ||a £ (t)||^2 
ll« e (*)lli» = KIU»> 



Z\ (age*?/*) - Six 



L- 



< ||a? - a e (t)|| L2 + |L £ (t) (e^ /£ - e **"W/ e ) 



<|L4 £ (i)|U 2 + ||a £ (i)|| 



L 2 



2s 



^Ct 2 + £\\*%t)\\ L ~ <-, 
2e e 



where we have used Theorem 15^4 



□ 



Corollary 5.6 (Local error for quadratic observables) . Let s > d/2 + 1 and \i > 0. 

Let <j) £ G L°°,a £ G H s with 

||^||l« HV^Uff.+i «S n, ||a £ || ffs ^/!. 

There exist C, Co > independent of e G (0,1] smc/i i/iai /or ^5 £ ^ Co, and 



\ z l u l\ 



\Slu e 



e"0| 

Im (eZlu e \7Ziu £ Q 



L 1 (R d )nL=°(R d ) 



- Im ( e5*u§V5^«g 



L 1 (R d )nL°°(R d 

Proof. Resuming the notations from the proof of Corollary 15 .51 we have 



\zy \ 



\sy \ 



,e|2 



K(<)l : 



and Cauchy-Schwarz inequality yields 



IKI 2 - l« £ W| 2 || L1 < ll«f - « £ WIIl 2 (IKIIl 2 + ||a £ (t)|M • 



The first part of the corollary then stems from Theorem 15.41 Similarly, 
Im (eSjSgV^ug) - Im (e^V^^) = |a?| 2 V^ E - \a £ (t)| 2 V <ff{t) 



+ slm(a £ t Va £ t ) - e Im [a £ (t)V a £ (t) ) . 

The second part of the corollary then follows easily from Holder inequality and 
Theorem 15.41 □ 
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6. End of the proof of Theorem 11.51 
6.1. Lady Windermere's fan. We denote 

(■?At\ n I 90 

{e } u 

To prove Theorem 11.51 we rephrase it in a more precise way: 

Proposition 6.1. Let s > d/2 + 2, 4> G L°°,a G H s , with V0 O G H s+1 , and T 
as in Theorem \1.5[ There exist v, 7, Afo, C\, Cq > such that for all e G (0,1], all 
^ At < At and all neN such that t n = nAt G [0, T], 

(6.1) I|V#||h. + |K||jj.-i 

(6.2) \\V<f> e n - V<f> £ (t n )\\ H s + ||4 - a e (t n )\\ Hs -i ^ 7 At, 

(6.3) ||V^|| ff3+1 + < e CiraAt ^ Co = e Cll/T , 

(6.4) IK-0 e (*n)||L- ^7At. 

Remark 6.2 (L°° bounds). The above result has an important technical conse- 
quence: the numerical solution it^ = a^e 1 ^™/ 6 is uniformly bounded in L°°(R rf ). 
In view of Proposition 13.21 the same holds for the exact solution u £ (t). Such in- 
formations are very delicate to obtain in general. Even in one dimension, the 
Gagliardo-Nircnberg inequality 

KIU- < V2\\u s \\% 2 \\d x u e \\% 2 

would not yield better than ||u e j|i=o < e -1 / 2 , because of the rapid oscillations 
present in u £ (4> e 7^ 0). Here, the uniform L°° estimates follow from the fact that 
a WKB regime is considered. 

Proof. The proof that we present follows essentially the lines of [HI Section 5]. 
Denote by (^ k ^j = (^f)' (tf^j numer i ca l solution, and 

k \ / ± 

g{n-k)At ( 9k 



,k 



«k 



In this proof, we omit the dependence of all the functions with respect to e. From 
Proposition 13. 2\ there exists R such that 

(tm) :=\\Vm\H^ + \Ht)\\H^R, Wg[0,T]. 

\ a \ L ); VH s+i xH s 

We prove Proposition 16. II by induction, with u = R + 5, <5 > so that the solution 
to p,12[) with data in the ball characterized by (|6.ip remains smooth up to time T 
(this is possible, since T < T max ). The estimates are obviously satisfied for n = 0. 
Let n ^ 1, and suppose that the induction assumption is true for ^ k ^ n — 1. 
We introduce the same telescopic series as in [T8], which is different from (ll.lOp . 
the latter being useful mostly when the problem (hence the splitting operator) is 
linear: 



^ ( 5 ) " ($£]J = E [S^-^Z^ (Jj - £(n-j-l)AtgAt ^ 



-1 



J 
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Noting the properties /„ = /" and f{t n ) = fn(f = ( l } or a ): we estimate 
\\V4> n -V(j)(t n )\\H° + \\a„ - a(t n )\\ H s-i 



n-l 

fe=0 
n-l 



(IIV0: 



fc+i 



,fc+l 



a n \\H>-^\ 



k=a 



g(n—k—l)At 
A/ 



rAt 



_ s (r 



-k-l)At 



S 



A/ 



VH s xH s 



For & ^ n - 2, 2£ 



V a fc+i 



and Proposition 13.31 yields, along with the 



induction assumption (all the norms are in VfP x H s ), 



?A/ 



s At 

< K(2R) 



-S 



A/ 



KjAt + R, 



<Kt k 
a(t k 

a(t k ) 



S 



A/ 



a(t k ) 



Htk+i) 

a{tk+i) 



which is bounded by R+S if < At ^ Ato <C 1. Up to replacing K with m&x(K, 1) 
we obtain that, for k ^ n — 1 and nAt ^ T, 



is controlled by 



fc-l)At 



rAt 



-fc-l)At 



s 



Ai 



VH°xH s 



K 



rAt 



-s 



A/ 



Using the local error estimate from Theorem 15. 4[ we infer, using 



g(n-k-l)At / g 



rAt 



_ g(n-k-l)At / g 



-At 



^ CK (At)' 



H- 



for some uniform constant C depending on Co- Therefore, 

||V</>„ - V(f>{t n )\\ H s +\\a n - a(t n )\\ H s-i ^ CTKAt, 

and we can take 7 = CTK, which is uniform in n and At, in order to get (|6.1[) 
and (|6.2p . Then (|6.3[) follows from Corollary 14.81 in view of (|6.ip and Sobolcv 
embedding, since we have assumed s > d/2 + 2. Finally, the L°°-estimatcs f)6.4[) for 
4>n are now straightforward (up to increasing 7), and are left out. □ 

Remark 6.3 (Nonlinear Schrodinger equation). We can now explain why Assump- 
tion [L2] is needed for the complete argument to work out. If we wanted to prove 
the analogue of Theorem II. 51 for, say, the defocusing cubic Schrodinger equation 

ied t u € + yAu e = \u £ \ 2 u e , 

then many results would still be available. In terms of the numerical scheme, the 
only change would affect the operator y l e : (|2.2[) would be replaced by 
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Working in H s for s > d/2, we see that unlike what happens under Assumption 1 1 . 2 [ 
4> e cannot be more regular than ag. On the other hand, the WKB formulation of 
the free Schrodingcr flow (|2.3[) induces a shift of regularity: if <fi £ is in H s for s 
large, then a £ must not be expected to be more regular than H s ~ 2 . Therefore, the 
splitting operator Z\ induces a loss of regularity, and this loss is iterated like T/ At 
times. It is this aspect which makes it hard to adapt Proposition 16. II to the case of 
the nonlinear Schrodingcr equation. 

6.2. Proof of Corollarv ll.6l Once Theorem ll.5l is available, we simply write, like 
in the proof of Corollary 15.51 



« 



o e (t n ))e < ^/ e + o e (t n )(e 



Taking the L 2 -norm, we infer 



-si 



L- 



< \K-a e {t n )\\ L 2 + \\a £ (t n )\\ L 2 



<t> e {t n ) 



and Corollary 11.61 is a direct consequence of Theorem 11.51 

6.3. Proof of Corollary 1 1.7t Corollary 1 1.71 also stems directly from Theorem ll.5 
by resuming the same computations as in the proof of Corollarv l5.6l 



7. Weakly nonlinear regime 

We now consider (jl.ip in the case a ^ 1, under Assumption 11.81 on the nonlin- 
earity. In view of the formal computations presented in Section [1.2| the case a = 1 
can be considered as the only interesting one, since no nonlinear effect is expected 
at leading order when a > 1. Since it is possible to treat both cases at once, we 
take advantage of this opportunity. 

The analysis in the case a ^ 1 being quite easier than in the case a = (even 
under Assumption 11.21 which is weaker than Assumption 1 1 .8[ ) . we shall simply un- 
derline the modifications to be made in order to prove Proposition ll.lOl bv following 
the same steps as in the proof of Theorem 11.51 

To characterize the exact flow in terms of WKB states, (|2.6|) is replaced by 
fa t e + i|V0 e | 2 = O; tf t=0 = <h, 

w 1 

[ d t a e + V<f • Va e + -a E A<j> £ = C-Aa e - ie a ~ Y ] (|a e | 2 ) a E ; af t=0 = a Q . 

Thanks to the assumption a ^ 1, the last term in the equation for a £ is not singular 
as e — > 0. More importantly, this is no longer a coupled system: the first equation 
is an eikonal equation, which we have analyzed in Section 14.21 

In the numerical scheme, the operator X*, corresponding to the free Schrodingcr 
flow, is the same as before, and analyzed in Section l4~2l On the other hand the oper- 
ator y* e can be modified, since the nonlinearity does not affect the rapid oscillations 
(as can be seen also from (|7.1|) ) . We recall that we now consider 

d t 4> £ = 0; cf>f t=0 = <p%, 

d t <f = -ie a - l f{\a°\ 2 )cf; af t=0 = a§. 
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We see that the possible loss of regularity pointed out in Remark 16.31 is not present 
here, since the regularity of <fi £ is not affected by the regularity of a £ . Also, working 
with a £ in H s for s > d/2 ensures that the analysis of Section B~2l can easily be 
adapted under Assumption 11.81 since a £ (t) = agexp (— is a ~ 1 tf(\a £ ) \ 2 )) . 

The main modification in the analysis concerns the local error estimate, since 
the statement of Lemma 15.21 must be revised. The operator A remains unchanged, 
and the operator B becomes 



D 

We compute successively 
B' 





-ie<*- l f(\a\ 2 )a 



n = -,>«-! ( ° 

b) \2fr (Re (ab)) a + h (\a\ 2 ) b + 2f> (\a\ 2 ) Re (ab) 



and 



where 



[A,B]r\ =is 





F(<f>,a) 



F(<j>, a) = V^> • V (/ (M 2 ) a) + -f (\a\ 2 ) aA0 - i|A (/ (\a\ 2 ) a) 

- a div /i (|a| 2 - ea div /i (Im (aVa)) 

- /2 (|a| 2 ) div (|a| 2 V0 + elm (aVa)) . 

The main point to notice is that if s > d/2 + 2, then F maps H s x H s to H s ~ 2 . 
Proposition 1 1 . 1 Ol then follows by resuming the same steps as in the proof of Propo- 
sition 16.11 



Appendix A. Linear Schrodinger equation 
Consider the (linear) Schrodinger equation with a potential, 

£ 2 

(A.l) ied t u £ + yAu e = Vu £ ; u e t=0 = itg, 

with V = V(t,x) <E R. We assume that V grows at most quadratically in space: 

Assumption A.l. V S L^ c ([0, oo) x H d ) is real-valued, and smooth with respect 
to the space variable: for (almost) all t ^ 0, x i— > V(t, x) is a C°° map. Moreover, 
it is at most quadratic in space: 

Va e N d , \a\ > 2, VT > 0, d^V £ L°°([0, T] x R d ). 

In addition, t h-> V(t, 0) belongs to L^ c ([0, oo)). 

Then for G L 2 (R d ), (|A.1|1 has a unique solution u £ g C([0, oo); L 2 (R d )), and 
its L 2 -norm is conserved, ||u E (i)||£2 = ||uo||l 2 for all t ^ 0; see e.g. [13] , We have 
more precisely: 

Proposition A. 2. Let k E N, V satisfying Assumvtion \A.1[ and Uq <G L 2 (R d ). 

Suppose in addition that Uq satisfies 

(A.2) HuSIIe* := sup (\\u £ \\ L2 + |||.x| fc ^|| L2 + || | £ V|^|| L2 ) < oo. 
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Then for all T > 0, the solution to (|A.1[) satisfies 

sup sup (\\u e (t)\\ L2 + \\\x\ k u e (t)\\ L * + |||eV|V(t)|| £a ) < oo. 

0<e^l te[Q,T] 

Proof. The key point is that the functions eWu £ and xu £ satisfy a closed system of 
estimates. Indeed, eV does not commute with the equation, and eVir satisfies 

iedt (eVu e ) + —A (sVu e ) = VeVu e + (eW) u £ . 

Similarly, 

f 2 



The standard L 2 estimate then yields 

\\eVu £ (t)\\ L 2 < ||eV^|| L2 + / ||(V^K(r)|| L2 dr, 

Jo 

||ot e (*)|| L 2 <; Hxugllia + / \\eVu e (T)\\ L 2dT. 
Jo 

Now under Assumption IA.I1 for T > fixed, we have the pointwise estimate 

\VV{t,x)u £ (t,x)\ ^ C(T) (1 + \x\) \u £ {t,x)\ , O^t^T. 
Recalling that the L 2 -norm of u £ is bounded, Gronwall lemma, applied to 

y{t) = \\e\7u £ {t)\\ L 2 + \\xu £ {t)\\ L 2, 
yields the proposition in the case k = 1. The general case follows by induction. □ 

Example A. 3. If Uq is of WKB type Q1.2p . or more generally (|2.1[) . with </>q a t most 
quadratic (in the sense of Assumption I A. ip . and oq £ H k D J r (H k ), then the above 
assumptions are fulfilled. Note however that Proposition IA.2I is valid for all time, 
and in particular after the formation of caustics, if any. 

Example A. 4. If 

«8(») = ^ (^) e ^ 9) - p/e > 

with q,p £ R d , 8 £ [0, 1], and ao £ S(R d ), then again, Proposition IA.2I is valid for 
all time. If 8 = 0, this datum is a particular WKB datum (with a linear phase). If 
8 = 1/2, this means that an initial coherent state is considered (see e.g. pT|). If 
8=1, the initial datum is concentrating at point q, which corresponds to a caustic 
reduced to one point (focal point; see [7]). 

Recall that if the splitting operators are defined by 

A = i-A; B = —-V, 
2 e 

then their Lie commutator is given by 

[A, B] = W- V + ^AV. 

With the norm ||u||e2 is defined in (|A.2|) . note the control 

||eW-Vu|| £a < \\u\\ B 2, 
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which follows from Assumption I A. ll By working with the norm rather than 

with the norm defined in Section ITTTl and used in [21 [H], the following result 

is a direct consequence of [12] and Proposition IA.2I 

Proposition A. 5. Let d 1, and V satisfying Assumvtion \A.i\ Suppose that 
||uo||s2 < oo. Then for all T > 0, there exist C, cq independent of e £ (0,1] such 
that for all At £ (0, cq], for all n £ N such that nAt £ [0, T], 



,A(\» , 



St 



where S^Uq = u e (t) in (|A.1|) . and Z* = e tB e 



tB„tA 



L 2 (n d ) e 
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